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Abstract 



The purpose of this master thesis is to simulate suitable nuclear reactions for a commissioning 
experiment, to be performed with the AGATA 7-ray tracking spectrometer. The main aim of 
the work is to find a reaction, which gives large Doppler effects of the emitted 7 rays, with as 
small contribution as possible due to the energy and angular spread of the nuclei emitting the 7 
rays. Inverse kinematics heavy-ion (HI) fusion reactions of the type (HI, 7), (HI, n) on proton and 
deuteron targets have been studied. Target effects were investigated using the program TRIM in 
order to determine the impact on the Doppler effects caused by energy and angular straggling in 
the target material. The cross sections of a large number of reactions of protons and deuterons on 
nuclei with mass numbers in the range A Ri 20 — 100 have been evaluated using the TALYS reaction 
code. The fusion-evaporation reactions, d (^^V, n)^^ Cr and d (^^Cl, n)^^ Ar were simulated in detail 
using the Monte Carlo code evapOR. The interactions in AGATA of the 7 rays emitted in these 
reactions were simulated using Geant4. The energy resolution of the 7 rays after 7-ray tracking 
and Doppler correction were determined as a function of the interaction position resolution of the 
germanium detectors. The conclusion of this work is that of the two reactions d (''^V, n)^^ Cr is 
more suitable for an AGATA commissioning experiment. 
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Chapter 1 

Introduction 



Within the field of nuclear physics one attempts to describe how the nucleus as a composition of 
nucleons can be held together, and how the nucleons interact both individually and collectively 
inside the nucleus. Nuclear models which describe the nuclear structure, can be tested and further 
developed by improving nuclear detection methods. The chart of nuclides, see figure [TTT] currently 
consists of about 3500 known nuclides. The green area of the figure [T]T] shows the nuclides located 
within the drip-lines. According to theoretical calculations less than half of the nuclides are known. 
In order to produce and study the nuclear structure of these unknown nuclides, great challenges 
face the development of precise detection instruments. This master thesis aims to design a commis- 
sioning experiment to be performed with the new 7-ray spectrometer. Advanced GAmma Tracking 
Array (AGATA) ^ . The experiment will be simulated and based on these calculations a possible 
reaction for the upcoming experiments in 2009 will be proposed. 
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Figure 1.1: Chart of nuclides. (Figure from Holifield Radioactive Ion Beam Facility, Oak Ridge 
http : // www . phy . ornl . gov I 



1.1 Nuclear spectroscopy 

Spectroscopy early became an effective way of revealing information about the atomic and sub- 
atomic quantal systems. In atomic physics, studies of emission spectra corresponding to the tran- 
sitions between electronic orbits, was a key element in understanding the atom. Nuclear physicists 
also use spectroscopic tools, by which one may retrieve valuable information about the inner struc- 
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ture of the nuclei. The energy levels (states) of the nuclei are denoted by their spin (J), and parity 
(tt), as (J'^). The de-excitation process of the states occurs either through emission of particles 
(e~,a ,P--.) or 7 rays. The 7 radiation emitted by the de-exciting nucleus spans an energy range 
from a few keV to tens of MeV. Gamma-ray spectroscopy plays a pivotal role in discovering new 
exotic nuclei far beyond the line of stability. Theoretical models can be evaluated and improved 
by testing them with the experimental data obtained by using 7-ray spectroscopic tools. Among 
the 7-ray spectroscopy instruments used so far, the high-purity germanium detector (HPGe) is the 
most successful one. By assembling a large number of such detectors in an array surrounding the 
position of the nuclear reaction, one can obtain efficient and precise detection of the emitted 7 
rays. 

One major factor contributing to the uncertainty in the measured 7-ray spectra is the Doppler 
effect. Considering in-beam experiments, the nuclei emitting 7 radiation are usually not at rest 
in the laboratory frame; hence the measured energy of the 7 rays is subject to an uncertainty 
caused by the Doppler effect. If not corrected for, this will lead to the broadening of the 7-ray 
spectra. The Doppler effect depends on two factors, the velocity of the nucleus and the angle 
between the direction of motion of the nucleus and the emitted 7 ray. Due to the second factor, 
the first interaction point of the 7 ray in the detector should be determined in order to minimize 
the Doppler effects on the broadening of the 7-ray peaks. One main advantage with new 7-ray 
spectrometers compared with previous instruments, is the 7-ray tracking technique, by which the 
first interaction point can be determined with much better precision than before, which enhances 
the Doppler correction capabilities. At present two such 7-ray spectrometers are being developed, 
AGATA 2^ in Europe and GRETA pi in USA. 

1.2 Motivation 

It is of great interest to study how one can compensate for the Doppler broadening of the 7-ray 
peaks in order to achieve a higher energy resolution with the AGATA spectrometer. The task of 
this work was to find a suitable nuclear reaction, in which the Doppler shifted 7 rays can be studied 
using one AGATA triple cluster detector in an experiment to be performed at INFN/LNL in Italy. 

One of the aims of this experiment is to use the Doppler effects to determine the interaction 
position resolution of the AGATA germanium detectors. Figure \L2\ shows a schematic view of the 
experimental setup. The reaction in study is 

A + a^B + b, (1.1) 

where A is the accelerated ion colliding with target a. The products after the collision are the 
particle b and the residual nucleus B, which emits the Doppler shifted 7 ray. 
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Figure 1.2: A schematic view of the experiment where the ion A from the accelerator coUides with 
the target a. 9b and 9b are the angle of emission of the residual nucleus and the emitted particle 
respectively, and 9^ is the angle between of the direction of the residual nucleus B and the emitted 
7 ray. 
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Chapter 2 



Theory 



2.1 Nuclear models 

There are essentially two classes of nuclear models. The first class is called independent particle 
models, where the Pauli principle restricts the collisions of the nucleons inside the nucleus, leading 
to a large mean free path. The shell model belongs to this class. The second class of models 
assumes strongly interacting nucleons inside the nucleus, resulting in a small mean free path. The 
nuclear liquid drop model belongs to this class. In general these models are called collective models 
as they attempt to describe phenomena that involve the nucleus as a whole, e.g. vibrations and 
rotations. 

In this chapter brief descriptions of the shell model and of the collective model are given. 
2.1.1 The shell model 

According to the Pauli principle the fermions cannot occupy the same quantum state. Therefore 
nucleons can orbit almost freely inside the nucleus, similarly to the electrons in the atom. In 
the atom, electrons are moving in a central potential caused by the positively charged nucleus. 
In the nucleus, the potential is created by the nucleons themselves through the nucleon-nucleon 
interaction. In the shell model, the nucleons move in a central potential usually parametrized as 

V{r)^—^^^ + VUr)t-^. (2.1) 

The first term in equation |2.1| is the Woods Saxon potential with the parameters: potential depth 
Vo, diffuseness parameter a, nuclear radius R — tqA"^/^ {A is the nucleon number and vq is the 
radius parameter). The second term is the spin-orbit interaction, where the spin {~s) and angular 
momentum ( £ ) vectors can be oriented either parallel or anti-parallel in the total nucleon spin 
quantum number j = i zL ^. 

The single-particle model is an extreme shell model form where the nucleons move freely in a 
central potential without interacting with each other. An odd-A nucleus in this model is composed 
of an inert even-even core, plus an unpaired nucleon. The last unpaired nucleon determines the 
nuclear properties, such as spin and parity. All known even-even nuclei have spin and parity 

— 0+ in the ground state. By introducing a residual force in the shell model, called pairing, 
one can explain the ground state spin and parity of even-even nuclei. The pairing interaction 
implies an attractive force between pairs of nucleons with opposite orbital angular momenta in 
identical orbits. The strongest attractive interaction occurs for nucleon pairs with largest possible 
quantum numbers rrij of opposite sign, i.e. ruj = +j and rrij = —j. The quantum number rrij is 
the projection of j on a quantization axis and can have values from —j to j. The spin and parity 
of each pair of nucleons in identical orbits with ±TOj is 0"''. The attractive pairing interaction leads 
to an increase of the binding energy for a completely paired nucleus with total spin and parity 0^ , 
compared with one with any number of non-paired nucleons. 

The shell model can explain the presence of magic proton (Z) and neutron (N) numbers in the 
nuclei. By including the spin-orbit term in the potential, the shell model gives the magic numbers 
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Z, N = 2, 8, 20, 28, 50, 82, 126, which corresponds to strongly bound nuclei at filled shells. 
Experimental results verify these theoretical shell model predictions. Doubly magic nuclei have 
magic proton and neutron numbers. The few nuclei which are doubly magic are strongly bound, 
like e.g. ''He. 



2.1.2 The collective model 

Two major modes of collective excitations are described within this model, namely vibrational and 
rotational excitations [3 . 



2.1.2.1 Vibrational excitation of spherical nuclei 

The nucleus experiences oscillations which could either change the size leaving the shape constant 
(breathing mode), or change the shape and leaving the density constant. The latter is more usual 
and can be explained by different modes of oscillations. The surface of a nucleus oscillating around 
a spherical equilibrium shape is written as 

(oo A \ 
1 + E E aA^Wi'AM(^:'^) , (2.2) 

where a\fj,{t) is the time dependent shape parameter, Rq is the radius of the sphere, and Yxf_i{9, (f) 
are the spherical harmonics. The multipolarity A determines the oscillation type, with A = 
corresponding to monopole oscillations. These oscillations have no angular dependency. Dipole 
oscillations occur when A = 1, causing the nucleus to vibrate around a fixed laboratory reference 
point. The value of A = 2 implies quadrupole oscillations in which the nucleus alternates between 
prolate and oblate shape. The energy levels of these vibrational modes are given by 

En,x = nLOxY^ (riA,, + 1/2) = hLUx (nx + ^ (2A + 1)] , (2.3) 

where N\ is the number of phonons in the nuclear vibration. A phonon is one quantum of vibra- 
tional energy, in analogy with the photon being the electromagnetic quantum. The frequency is 
given by uj\. An increase of the number of phonons implies higher energies, for a vibrational mode 
A. These energies will form a ladder with increasing hujx- This spectral phenomenon is called a 
vibrational band 4 . 



2.1.2.2 Permanently deformed nuclei and rotational excitation 

The spherical symmetry is broken when moving away from closed shells. Due to the pairing effect 
mentioned in 2.1.1| the nucleons outside closed shells, will occupy the largest possible ±TOj values, 
which leads to a non-spherical shape. The existence of permanent deformations will open up the 
possibility for nuclear rotations, which are not present in spherically symmetric nuclei. Due to 
the rotation the nucleus feels a centrifugal stretching since it is not a rigid body. In analogy with 
vibrational bands, rotations form rotational bands. 



2.2 Nuclear reactions 

Nuclear reactions can be elastic, if the target and the projectile are in the same states before and 
after the reaction (usually in their ground states). When excited states are formed or new particles 
are created, the collision is inelastic [5 . When nuclei collide with each other many reaction types 
can occur. Direct reactions and compound reactions are two important types of reactions. 

2.2.1 Direct reactions 

In a direct reaction the projectile can undergo stripping in which case some of the nucleons are 
transferred from the projectile to the target. An example of a stripping reaction is (d,p). A 
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knockout reaction can also occur, in which the incident nucleus knocks out nucleons from the 
target nucleus, as in a (p,np) reaction. Direct reactions occur on a very short time scale, 10^^^ to 
10~^^s 3 , and do not change the structure of the projectile and the target dramatically. 

2.2.2 Compound nucleus reactions 

According to |5j, compound nucleus reactions occur when two nuclei form an intermediate state. 
The compound nucleus formation between nuclei a and A is depicted as 

a + X* ^ R* + J2 particles, (2.4) 

where X* is the excited compound nucleus and R is the excited residual nucleus. The timescales of 
these reactions are relatively long, roughly 10^^^ to 10~"'^^s. The information about the constituents 
before the collision is permanently destroyed, hence the output decay channels is independent of 
the compound nucleus formation. If the same compound nucleus is reached through another 
reaction, say c + C —^ X*with different colliding nuclei, the probability of the decay channel 
X* —rR* + '^ particles remains the same. This hypothesis is called the independence hypothesis 
and it implies that it is possible to separate the cross section a, into the formation cross section 
a^^ for projectile a and target A, and the branching ratio Gg^ (CN=compound nucleus) for the 
decay to particle b and nucleus B: 

a = a^^iE)Gg^iE). (2.5) 

The cross section term tT^J^(i?) is the probability of the compound nucleus formation at the center 
of mass energy E. As the compound nucleus must decay, the sum of all the partial decays is 
^Gp^(i?) = 1. The excitation energy E^-^ of the compovmd nucleus can be expressed as 

i 

= y , , AT -^'^b + [ma + TTlA - mcN] (2.6) 

where i?iab is the kinetic energy of the projectile in the laboratory frame and the masses m^, rriA 
and nicN are the atomic masses of the projectile, target, and compound nucleus, respectively. The 
compound nucleus de-excites through particle evaporation. 

2.2.2.1 Compound nucleus formation 

Depending on the collision type and the impact parameter, high orbital angular momenta may 
be transferred to populate the compound nucleus. The maximum amount of angular momentum 
transfer in the compound nucleus reaction, is denoted with £cr [S]- By considering other reaction 



types, higher £ values can be reached. Figure 2.1 shows the differential cross section da/dl as a 
function of £. 
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Figure 2.1: Differential cross section as a function of orbital angular momentum. 



Direct reactions occur for peripheral collisions implying high angular momentum transfers. 
When the total angular momentum transfer £gr is reached, the cross section drops. The sharp 
cut-off model states that this drops instantly when £g^ is reached. This assumes a well-defined 
nuclear radius and a behavior similar to the classical rigid-body collision. But since experiments 
show that the nuclear radius is diffuse, and the border is not well-defined, the collision may not 
fall drastically when is reached. Instead of this sharp cut-off model, a smooth drop-off model 



is showed in figure 2.1 [5", 



The angular momentum transfer is essential in order to understand the probabilities for nuclear 
reactions. According to the statistical model, for heavy-ion collisions, the reaction cross section for 
the compound formation depends on the transfer of orbital angular momentum £ as [B] 



■Kk' 







{2i+l)TfPt 



CN 



(2.7) 



where is the transmission coefficient, k is the wave number, and Pf'^ is the probability of 
the ^-wave to enter the compound nucleus. In compound reactions the transmission coefficient 
Ti approximately equals unity for relatively high angular momenta. The cross section can be 
simplified into 



CN 



(2.8) 



The sum runs from orbital angular momentum until the critical value 
compound nucleus formation cross section rapidly drops to zero. 

The cross section for final formation and decay is finally given by [7] 



in which value the 



Jc ^ j\ cj y s=\i-s\e=\j^-s\ ^ ' ) 

(2.9) 

where (JaA {E) is the compound nucleus formation cross section, dependent on the entrance channel- 
energy E. The decay width Vi is for a specific final state i in the exit channel. The target spin is 
given by /, and S is the channel spin. The projectile spin is represented by s. 
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2.2.2.2 Compound nucleus decay 



When the compound nucleus is formed it populates energy levels high above the ground state. 
The nucleus then tries to reach the ground state through an evaporation process in which mainly 
light particles (neutrons, protons, a particles) are subsequently emitted in a so called cascade, see 
figure [2T2I While charged particles, such as protons and a particles, have to overcome the Coulomb 
barrier in order to escape the excited nucleus, neutral particles do not. This will be reflected in 
the energy distribution of the emitted particles: neutrons are emitted with lower kinetic energies. 



See figures 6.14 and 6.21 for a comparison of energy distributions of neutrons and alpha particles. 




Figure 2.2: In the left figure a population scheme is shown, in which compound nuclei are created 
with a given excitation energy E* . Particle evaporation releases most of the energies, followed by 
7 ray cascades to finally reach the ground state. In the right figure the angular distribution of 
neutrons and a particles are schematically presented, indicating an increased anisotropy for the a 
particles. 



Alpha particles carry away larger values of angular momentum from the compound nucleus 



than neutrons and protons, as shown in figure 2.2 This follows since 



£ = ~r*x'p* = 7^x ml? (2.10) 
and both the mass to, and the velocity v, is higher for a particles than for neutrons and protons. 



Energy distribution of evaporated particles 

The energy distributions of the evaporated particles have a Maxwell-Boltzmann shape, according 
to the spin dependent statistical model. Following the derivation of Lang and Lecouteur 8 , the 
probability of emission is 



Wp{E*,J,E*^,Jj^) 



2S, 



p ' 



p{E*,J) ■ 



(2.11) 



where the emitted particle has a reduced mass /i, energy E-p, and spin 5p. The inverse compound 
nucleus production cross section is CTj(j^(i?*), and p{E* , J) and p{E^, Jr) are the level densities of 
the initial and residual nuclei at spin and excitation energy E* , J and E'^, Jn, respectively. 

Angular distribution of evaporated particles 

Based on the work of Ericson and Strutinsky ^ , the orbital angular momentum of the incident 
particle is orthogonal to the direction of the beam. The total angular momentum aligns preferably 
in the plane perpendicular to the beam axis. The angular distribution is symmetric with respect 
to the angle 90° relative to the direction of the incoming beam. According to [6j[9], the angular 
distribution of the evaporated particles can be written as 



W. 



477 



(2fc)! 

(2'^fc!; 



J2k 



iJi 

„2 



P2fe(cOS 6'p). 



(2.12) 
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The spherical Bessel functions j2k are of order 2k and 9p is the angle between the incident beam 
and the direction of the emitted particle in the center of mass system. The spin cut-off parameter of 
the residual nucleus is given by cTf . The Legendre polynomial P2k is even, thus there is a symmetry 
around 90°. Figure 2.2 illustrates that a particles are emitted mostly at large (~ 180°) and small 
0°) angles. The reason for this is that they remove large values of the angular momentum £ 
when evaporated from the compound nucleus. For neutrons the evaporation is nearly isotropic. 
See figure |6.12| for an angular distribution of neutrons obtained by the simulations performed in 
this work. 
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Chapter 3 

Nuclear spectroscopy methods 



3.1 Interaction of 7 rays with matter 

When 7 rays enter matter, three interaction types can occur, photoelectric effect, Compton scat- 
tering and pair production. 



3.1.1 Photoelectric effect 

When irradiating a metal surface with light above a certain frequency, an electric current will flow. 
Einstein found that electrons are removed from the metal and that the photon is absorbed by the 
atom [TD]. The photo-electron can only be ejected when the photon energy is higher than the 
binding energy of the electron. The kinetic energy for the photoelectron is given by 

' = i?^ - (3.1) 

where is the photon energy and Be is the electronic binding energy. The photoelectric effect 
dominates at low (see figure 3.11. 



3.1.2 Compton scattering 

When the photon energy increases, the probability of photon scattering with electrons increases (see 



figure 3.11. According to Compton in his article from 1923 |TT], the observed photon wavelength 
A increases after the interaction, leading to a decreasing energy, since E = hc/X. Recoil energy 
is transferred from the photon to the electron. Due to energy and momentum conservation, one 
arrives at the Compton scattering formula which gives the energy of the scattered 7 ray 

^ l + (S^/mc2)(l-cos0c)' ■ 

Here E^ is the initial energy of the 7 ray, and m is the mass of the electron, the scattering angle 
of the photon is 9q. 



3.1.3 Pair production 

If the photon energy is large enough, pair production may occur, i.e. creation of an electron- 
positron pair, as 7 ^ e^ -f e~. The minimum energy required is 1.022 MeV. Due to momentum 
and energy conservation this process is only possible if it occurs nearby a particle with non-zero 
rest mass, e.g. a nucleus [12 . 
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Figure 3.1: Domains of electromagnetic interaction with matter. [T^] 



3.1.4 Attenuation coefficients 

The intensity of a beam passing through matter, decreases exponentially 



(3.3) 



where Jq is the incident beam intensity, /Lt is the total attenuation coefficient and 5 is the material 
thickness. The probability for a photon removal per unit length is given by /i. The total attenuation 
coefficient is the sum of the attenuation coefficients for the photoelectric effect (r) , Compton 
scattering (c) and pair production (k) 



/i = T + c + K. 

The inverse of the attenuation coefficient is equal to the mean free path of the 7 rays 

1 



An 



(3.4) 



(3.5) 



3.2 Interaction of ions with matter 

Transport of ions in matter is mainly affected by electromagnetic interactions with the electrons 
of the material. Rutherford scattering, i.e. scattering between nuclei, is rare in comparison with 
the interaction with the electrons. The energy loss of the ion per distance traveled is given by [13] 



dE 
dx 



47rer 



mv^A 



In 



2mv 



iii-{v/cy 



(3.6) 



where A, Z, z, E and v are the mass number, atomic number, electric charge, energy and velocity 
of the ion respectively, m is the electron mass, p the density of the target, / is the mean excitation 
potential of the target, and A'pAvogadro's number. 



3.3 Doppler Effects 

Doppler Effects are an important concept for nuclear spectroscopy. Gamma radiation emitted by 
a nucleus in motion, appears to have a different energy than the actual transition energy. The 
frequency of the 7 rays is lower if the nucleus moves away from the observer, and higher when it 
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moves towards the observer. Figure [3T2| demonstrates the change of frequency /, which is observed 
by the receiver as /'. 




Figure 3.2: Illustration of the Doppler shift. 



The relativistic version of the Doppler shift formula for a source moving with velocity v is jl4j 



1 ^ ^ 

^^^^-^n-^cos^,' (^-^^ 

where Ej and E-y^ (E — hf) are the measured 7-ray energy and the 7-ray transition energy, 
respectively, and 9 is the angle between the receiver and the direction of motion of the source. The 
velocities in the simulations of the AGATA commissioning experiments are below 10 % of the speed 
of light, hence the relativistic factor -y/l — v'^/c^ is very close to unity. In this case a simplified 
non- relativistic approximation of equation |3.8| is appropriate to use: 



E — E 

^1 — I^jo ^ ^ 



(1 + ^cos6'.^) . (3.8) 



3.4 Gamma-ray spectroscopy 

When photons enter a 7-ray detector they may deposit all or part of their energy in the detector. 



Figure 3.3 shows how different interactions contribute to a typical 7-ray spectrum. 



PliDtopaaX 




Figure 3.3: Three 7-rays enter the detector in the left figure, a) Two Compton scatterings followed 
by a photoelectric effect, b) A pair production followed by single escape, c) Pair production 
followed by double escape. The figure on the right shows the main characteristics of a 7-ray 
spectrum. 
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Photons depositing all their energy in the detector give rise to the full-energy peak (photopeak) 
while the ones that do not deposit all their energy, but escape after one or more Compton scat- 
terings, contribute to the Compton continuum. The scattering angle in equation 3.2 determines 
the energy deposited. Back-scattering of the photons forms the Compton edge in figure 3.3 



amount of 7 rays that escape from the detector depends on the size of the detector compared to 
the mean free path of the photons in the detector. 



3.4.1 HPGe detector and energy resolution 

The high-purity germanium (HPGe) detector is undoubtedly the most successful instrument for 
spectroscopic studies of the nuclear structure. The superior performance of the HPGe detector is 
mainly due to its excellent energy resolution, which allows the distinction of 7 rays with close lying 
energies. The energy resolution is characterized by the full width at half maximum (FWHM) of 
the 7-ray peaks. Two factors contribute to the FWHM, here labeled W [15] : 

W= ^Jwf + W^. (3.9) 

The intrinsic width Wi is caused by the properties of the detector itself and consists of three terms: 
the electronic noise Wen, which is independent of the measured 7-ray energy, the statistical factor 
Wsf, which depends on the square root of the energy, and the incomplete charge collection Wcc, 
which is proportional to the energy, 

= ^/WL + W^^+Wl- (3.10) 

The second factor Wd is due to the Doppler effects and it is non-zero only if the 7 ray is emitted 
by a moving source. The Doppler factor consists of three terms: 

Wd = ^Wl+Wl + Wl. (3.11) 

Wq-^ is due to the broadening caused by the uncertainty in the determination of the emission angle 
of the 7 ray, 9^. This term depends on the position resolution of the first interaction point in the 
detector. The terms Wg^ and W^^ are due to the uncertainty in the determination of the angle 
and velocity, respectively, of the residual nucleus. 



3.5 AGATA 

AGATA is a 7-ray spectrometer which in its final A-k configuration will consist of 180 HPGe crystals 



assembled in 60 triple cluster detectors (see figure 3.4). Each HPGe crystal is of type closed-end 
coaxial. The crystals are 9 cm long, have circular shape with a diameter of 8 cm at the rear end, 
and a hexagonal shape in the front. The crystals are also electrically segmented into 6 longitudinal 
and 6 azimuthal segments which together with the core gives 37 separate signals per crystal. The 



crystals are encapsulated in aluminum cans (see figure 3.4 1 




Figure 3.4: The AGATA HPGe detectors [2]. 
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AGATA will use advanced digital electronics and pulse-shape analysis (PSA) techniques to 
determine the 3D position and the energy deposition in each 7-ray interaction point. An example 
of the PSA technique is shown in figure 3.5 The interaction points determined by the PSA are then 
fed into a 7-ray tracking algorithm, which is based on the Compton scattering formula, equation 
|3.2| and by which the full-energy and the first interaction point of each 7 ray is reconstructed. 
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Figure 3.5: Illustration of the pulse shape technique. In this event a hit in segment 4 can be seen. 
A signal is registered from the core and two mirror charges are seen in segments 3 and 5 [IB] . 



The improvement to be expected by AGATA compared to earlier Ge-detector arrays, not based 
on 7-ray tracking, will result in a larger full-energy efficiency and much better Doppler correction 
capabilities. With previous detector arrays one could only determine which detector was hit by a 
7 ray. Thus, one had an opening angle corresponding to the diameter of the detector, typically 
~ 5 cm. In AGATA the position resolution will be about 5 mm which will lead to a much smaller 
value of Wa 



in equation 



3.11 



and, thus, too much narrower 7-ray peaks. 
AGATA will serve as a powerful tool to understand nuclear structure far from the /3-stability 
line. Several detection difficulties such as Doppler broadening will be possible to compensate for 
using the AGATA spectrometer. By having a better resolution of the full-energy peaks, weak 
transitions in exotic nuclei can be studied [2]. Figure 3.6 displays the predicted performance of the 
AGATA demonstrator compared to the CLARA spectrometer. 
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Figure 3.6: The predicted performance of the AG ATA demonstrator (15 HPGe crystals) compared 
to the CLARA spectrometer, which consists of about 100 non-segmented HPGe crystals. The 
FWHM at a 7-ray energy of 1 MeV is more than a factor of two better with AG ATA (2.7 keV and 
2.9 keV) at a target to detector front face distance of 24 and 14 cm, respectively, than with the 
full CLARA array (6.2 keV) [17]. 
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Chapter 4 

Requirements 



The selection of reactions to be studied in this work was based on the following list of criteria. 

• Large Doppler shifts: This is required in order to study the interaction position resolution of 
the AGATA detectors. The velocity of the residual nucleus should be higher than 5 % of the 
speed of light. 

• Large reaction rate: The reaction rate is given by 

R--^^ (4-1) 

where / is the ion beam current, q the electric charge of the ions, t and p the target thickness, 
and density, respectively, a the cross section, ttia the atomic weight, and A^a Avogadro's 
number. 

• Available ion beams at LNL: The available beams at the LNL accelerator laboratory are 
listed in appendix [B| 

• Competing decay channels: Different particle evaporation channels compete. Proton and a 
particle evaporation give large recoil energies to the residual nucleus, since they need high 
kinetic energy to penetrate the Coulomb barrier. Neutrons do not face the same situation, 
hence the average recoil energy given to the nucleus is smaller. 

• Position of 7-ray emission: In order not to have 7-ray emission inside the target, the life 
times must be longer than the time it takes for the ions to pass the target material. To have 
a negligible influence on the Doppler effects, the life times must be short enough so that the 
decay occurs within 1 mm after the target. 

• Energy and angular straggling of the ions in the target should be as small as possible. 

• Gamma-ray energies: A minimum 7-ray energy of 500 keV is required. Lower energies will 
give smaller Doppler shifts, which are harder to study. 

• Strong 7-ray transitions: In even-even nuclei a strong 2+ 0^ transition is often observed, 
therefore even-even final nuclei are preferred. 

In this work two reaction types were studied in order to find a suitable reaction: a) proton capture 
p (A", 7) or deuteron capture d (A, 7) and b) fusion evaporation with emission of one neutron 
followed by 7 rays, for example p(A, n7) or d(A, n7). Both type of reactions were of inverse 
kinematics type, with a heavy ion beam on a proton or deuteron target. 
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Chapter 5 

Simulations 



5.1 TALYS 

The code used to evaluate the cross sections was TALYS, a comprehensive nuclear reaction modeling 
code [TS]. The projectiles and ejectiles in TALYS can be of type 7 rays, neutrons, -'^"^H, '^~'*He, 
with energies in the range 0.1 to 200 MeV. The mass number of the target can be in the range A— 12 
to A=339. Note that by TALYS it is not possible to simulate the reactions in inverse kinematics. 
TALYS was, however, used for calculating the cross sections of the reactions, which are independent 
of the reaction being of normal or inverse kinematics type. In this work the results obtained with 
TALYS were in several cases compared to experimental data taken from the National Nuclear Data 
Center. 

The total probability of a reaction to happen, i.e. the total cross section dtotj is the sum of the 
elastic CTc and the inelastic cross section (Tjc: 

C^tot = O-o + (Tic. (5.1) 

The inelastic cross section is a sum of the compound nucleus, pre-equilibrium and direct cross 
sections. The compound nucleus formation cross section a^, is defined as the difference between 
the reaction CTioac and the direct cross section (Tdiroct: 

Ccf = frcac fdircct- (5-2) 

The output of TALYS contains the calculated cross sections for non-elastic direct reactions, non- 
elastic and elastic compound nucleus reactions, and non-elastic pre-equilibrium reactions. 

5.2 Target effects 

The target that will be used in the experiment has an effect on the velocity and direction of the 
ions. The program used for studies of the target effects was TRIM [19 , which is a Monte Carlo 
based program that simulates collisions between incident ions and atoms in the target material. 
Angular and energy distributions of the incident ions are obtained from the TRIM simulations. 
The target material can be chosen from a list of chemical substances and compounds. 

5.2.1 Sources of target effects 

When the ions enter the target material they are deflected due to the electromagnetic interaction 
primarily with the electrons. Thus, the ions will loose energy and deviate from the beam direction, 
which leads to so called energy and angular straggling of the ions. There are two main aspects of 
the target effects in the simulations. 

1. Position of the reaction in the target. If the reaction occurs in the beginning of the 

target, the residual nucleus will experience energy and angular straggling, whereas if the 
reaction occurs at the end, the incident ion instead of the residual nucleus will lose kinetic 
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energy. The difference both in fcinetic energy and atomic number of the incident ion and 
residual nucleus leads to the difference in straggling. 

2. Neutron energy and angular distribution. After the reaction the excited compound 
nucleus is formed, and particles are emitted. Depending on the energy and angular distribu- 
tion of the particles, the residual nucleus will obtain different kinetic energies. The recoil due 
to the neutron evaporation will either boost the nucleus or slow it down. Since this happens 
inside the target material the energy and angular straggling will change after the reaction. 
Two extreme cases were investigated, one where the neutron is emitted at about 0° relative 
to the ion beam, and the other where it is emitted in at about 180°. 



5.3 Fusion-evaporation simulation 

The fusion-evaporation simulations were done using the Monte Carlo code evapOR [2D]. This 
program allows for fusion of two arbitrary nuclei, followed by evaporation of particles which can 
be of 7 different types (n, p, d, a, t,"^ He,^ Li) as well as emission of 7 rays. It uses the statistical 
model for compound nucleus reactions, described in section [2. 2. 2 [ 

In evapOR the velocity of the emitted particle in the laboratory system, l^B' obtained by 
adding the velocity of the compound nucleus in the laboratory system v q-^ and the velocity 
of the particle in the center of mass system Qjyj: 

" p,LAB CNXAB + "p,CM- (5-3) 

The transformation of the particle emission angle between the center of mass m and laboratory 
systems is given by 



tan 9 



sm ( 



p,LAB 



'p,CM 



(5.4) 



1/2 



(5.5) 



cos f p,CM 

where 

' mamb 

mATriB E + ( 

Here Q denotes the Q- value of the decay, ma, rnx, rn^,, me are the masses of initial and final 
products, and E the center of mass energy of the impinging particle. 

An output file with the format shown in table |5.1| is produced by evapOR. 





Index 


Energy 


vx 


vy 


vz 


Start new event 


$ 










Compound nucleus (Z,N) 


101 (18,39) 


63090 








1 


Emitted particle (q) 


7 


3615.4 


-0.81693864 


-0.13911782 


0.55969405 


7 


1 


2169.5 





1 






Table 5.1: Example of the .aga output file produced by evapOR. In the output shown here, the 7 
ray is emitted in the y direction. The parameters vx, vy, vz are the normalized (| v |= 1) velocities 
in Cartesian coordinates in the laboratory system. The indices 101, 7 and 1 denote the compound 
nucleus, an alpha particle and a 7 ray respectively. Energy is the kinetic energy of the particle in 
the laboratory system or the Doppler shifted 7-ray energy, all given in keV. 



The event shown in table |5.1| starts with the compound nucleus. The listed 7 ray was not 
generated by evapOR but was added artificially (21j. The energy of the 7 ray is Doppler shifted 
according to the non-relativistic Doppler formula of eq. (3.7). 



5.4 Simulation of 7-ray detection 

The complex geometry of the AGATA detectors have been implemented in the Geant4 simulation 
package [251 US]- In this package, one can choose between using a single HPGe crystal or any 
number of triple clusters up to the full Att AGATA array. 
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One can either generate particles in Geant4 or get them from an external event file. In the 
present simulations event files from evapOR was used as input to Geant4. The output of an 



ordinary run using an input event file similar to the output shown in table 5.1 is shown in table 
K2\ 







niimi^y 1 v/c 




y 




Start new event 


-100 










Compound nucleus velocity 


-101 


v/c=0.05oo9 


ON 


ON 


IN 


Position of CN 


-102 













Emitted particle (a) 


-7 


3303.4 


-0.04844N 


-0.89261N 


-0.44821N 


Residual nucleus 


-101 


v/c=0.06147 


0.00175 N 


0.03233 N 


0.99948 N 


Position of RN 


-102 













Emitted particle 7 


-1 


2172.300 


ON 


IN 


ON 


Interaction point 1 (detector 2) 


2 


383.495 


-0.373 


43.248 


-0.426 


Interaction point 2 (detector 2) 


2 


32.573 


-1.477 


43.844 


-1.851 


Interaction point 3 (detector 2) 


2 


34.293 


-7.208 


49.511 


-2.958 


RN after emitted 7 


-101 


v/c=0. 06147 


0.00175 N 


0.03133 N 


0.99951 N 


Position of RN 


-102 













RN after emitted 7 


-8 


66971.702 


0.00175 N 


0.03133 N 


0.99951 N 



Table 5.2: Example of an output file of the AGATA Geant4 package. N denotes that the velocity 
of the residual nucleus is normalized to 1. Energies are given in keV and positions in cm. RN 
stands for the residual nucleus. 



The events start (Index=-100) with the compound nucleus velocity divided by the speed of light 
in vacuum, c, and its direction (-101), followed by the position of the CN (-102), the evaporated 
particle (-7), and the residual nucleus (-101, -102). The nucleus then emits the Doppler shifted 
7-ray (-1), which interacts three times in detector number 2 (Index=2). After this the v/c and the 
direction (-101) and position (-102) of the residual nucleus is listed. Finally the kinetic energy and 
the direction of the residual nucleus is given (-8). 



5.5 Gamma-ray tracking 

The next step is to track the 7 rays [TB]. In this work the MGT tracking program was used. 
This tracking code uses the Compton scattering formula to determine the path of the 7 ray in the 
detectors and gives the full energy and the first interaction point of the incoming 7 ray. The first 
interaction point provides the angle 9^ needed for the Doppler correction according to equation 

In the planned AGATA commissioning experiment, no ancillary detectors for the detection 
of the residual nuclei will be used, therefore the velocity vectors of the residual nucleus will be 
unknown experimentally. Instead all the residual nuclei are assumed to have the average velocity 
and direction along the beam axis. 
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Chapter 6 

Results and discussions 



6.1 TALYS 

In table 6.1 the studied reactions are summarized. Following the requirements in chapter |4] most 
of the ion beams at LNL [appendix |b] were excluded from further study. In most reactions, 
competition between different reaction channels were found. The 7-ray energy was also in many 
cases lower than required. The proton and deuteron capture reactions had too small cross sections 
to be useful for the planned AG ATA commissioning experiment, hence only fusion evaporation 
reactions were studied. 



p - Reactions 


Result 


d - Reactions 


Result 


19F4 


P 


a channel dominates 


27A1 + d 


Competing p & n channels 


63Cu 


+ P 


p channel dominates 


63Cu + d 


Competing n, p & channels 




f P 


n ch. dominates, odd-even nucleus 


Ag + d 


Channels competing 


27A1- 


f P 


a &i p channels competing 


28Si + d 


p channel dominates 


31P4 


P 


p channel dominates 


37C1 + d 


n channel dominates 


35C1- 


fP 


p channel dominates 


siy + d 


n channel dominates 


58Ni. 


f P 


p channel dominates 


i™Mo + d 


p channel dominates 




f P 


p channel dominates 


94Mo + d 


p channel dominates 




-P 


p channel dominates 


«iBr + d 


p & n channels competing 


+ 


P 


p channel dominates 






33s + 


P 


p channel dominates 






65 Cu 


+ P 


n ch. dominates, bad ^ energy 







Table 6.1: Proton and deuteron reactions studied by TALYS. 
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Figure 6.1: Example of two (p,n) reactions, ^^Fe(p,X) and ^-^P(p, X). {a,— a particle, g—j ray). 
The energy of the incident proton is given in the laboratory system. 



In the left plot of figure 6.1 the °°Fe(p,p) Fe reaction channel is dominant. In the right plot 



of figure |6.1| the a emission dominates up to ~ 3 MeV at which energy proton emission starts to 
become important. 



Two examples of deuteron induced reactions are presented in figure 6.2 To the left, ^Si(d, X) 



where the proton channel dominates. In figure [6^ ^^Al(d, X) is shown, where there is a competition 
between the neutron and proton channels . 



600 



500 



400 



300 



200 



100 




300 



250 



200 



150 



100 



50 



Cross section for 27Al(d, X) 

27Al(d,d)27Al 

27Al(d,p)28Al 

27Al(d,a)24Mg 

27A1 (d, t) 25Mg 

27Al(d,n)28Si 

27Al(d,g)29S 



2 3 
E,[MeV] 




2 3 
Ed [MeV] 



Figure 6.2: Example of two (d,n) reactions, ^^Si(d,X) and ^^Al(d,X). The energy of the incident 
deuteron is given in the laboratory system. 



Results obtained from the TALYS simulations were verified against experimental data. Three 
examples of such comparisons are presented in figures |6.3| and |6.4[ 
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Figure 6.3: In the left figure the reaction ^^Cu(d,n) is compared to measurements reported in ref. 
[131111123111]. In the right figure 5'^Fe(d,n) is compared to experimental results from ref. [2^150] . 
The incident proton and deuteron energies are given in the laboratory system. 
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Figure 6.4: The ^^V(p,n) reaction compared to experimental results from ref. [5T1 [ 
of the incident proton is given in the laboratory system. 



The energy 



For the ^'^Cu(p,7) reaction (figure 6.3 left), four independent measurements were plotted to- 



gether with the calculation of TALYS. The differences are small at low energies but reach a factor 
of 2 around the peak at iJp « 4 MeV. In the plot on the right hand side of figure 6.3 the experimen- 



tal values are close to the results of the TALYS calculation for the ^"'Fe (d, n) reaction for energies 
up to about 14 MeV. At higher energies the measurements of ref. [15] differ from the TALYS 

the reaction V (p, X) is simulated and compared with measurements 



6.4 



calculations. In figure 

from ref. |311 132j . Both the neutron and 7-ray emission channels are close to the results obtained 
from measurements. 



6.1.1 Selected reactions 



In table 6.1 two reactions were shown with a dominating neutron emission channel. These are the 
reactions that were further studied. Both reactions leads to an even-even residual nucleus, and 
they have appropriate life times. The first reaction, '^^Cl (d, n) ^^Ar, leads to ■^^ Ar which has a first 
excited state 2+ at 2168 kcV. The high 7-ray energy which decays to the 0+ ground state is very 
good for the Doppler-effect measurements. The cross-section calculation from TALYS is shown in 
the upper plot of figure 6.5 The reaction ^^Cl(d, n)'^*Ar is clearly dominating. 
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The second reaction, ^^V (d, n) ^^Cr, is for the same reason a suitable reaction. The first excited 
state has an excitation energy of 1434 keV. The calculated cross section for this reaction is shown 
in the lower plot of figure [63) 



400 M I I — I I I I — II — I I I — I I I I — I I I — II — I I I — I I I I I — I I — I I 




EJMeV] 



Figure 6.5: TALYS cross-section calculations for •^^Cl(d,n)Ar and '^^V(d,n)Cr. The incident 
deuteron energies are given in the laboratory system. 



The cross section of the reaction ■^''Cl (d, n)'^^Ar reaches a maximum of about 350 mb at a 
deuteron energy of about 4 MeV. The cross section of the ^^V(d, n)^^Cr reaction has a local 
maximum at about Ej^lab = 3.4 MeV and it increases at higher energies without competition 
with other channels. Based on the TALYS results, both reactions, d -l-'^^ CI and d V, are 
suitable for further studies. The energies chosen are Ed = 3.8MeV and = 3.2MeV for the 
d+'^^Cl and d -1-^^ V reactions respectively. In inverse kinematics reactions thus correspond to 
£■37^^ = 70.0 MeV and £^51^ = 81.6 MeV, respectively. 



6.1.2 Comparison between evapOR and TALYS 

To verify the calculations, a comparison between TALYS and evapOR was made. The results for 



the d + V reaction are presented in figure 6.6 



25 




Figure 6.6: Comparison of cross sections obtained by TALYS (lines), and evapOR (symbols) for 
the d V reaction. The incident deuteron energy is given in the laboratory system. 



As seen in figure 6.6 large discrepancies were obtained. Below E « 3.4 MeV the TALYS and 
evapOR results are closer. Above 3.4 MeV the cross sections are similar in particular for the 
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V (d, n)^^ Cr reaction. 



6.2 Studies of target effects 

In the proposed commissioning experiment, an inverse kinematics reaction will be used. This 
means that for the two selected reactions, d CI and d V, deuterium would be the target 
nucleus. One possible target is to use a thin plastic foil (e.g. polyethylene) enriched in deuterium. 
Such foils do, however, only allow very weak beam currents without the foil getting destroyed 
(melting) . Instead it is proposed to use a deuterated Ti foil which can contain about 1 deuterium 
atom per titanium atom [55]. The chosen target thickness is 200 /im/cm . No deuterated titanium 
compound is available in TRIM therefore pure titanium was chosen as target in the simulation. 
The difference in energy straggling due to this is expected to be negigible. No nuclear reactions 
will occur between the '^^Cl or ^-'^V beams and the Ti nuclei of the target (the energies are well 
below the Coloumb barrier). 



6.2.1 Energy straggling 

The target studies were based on the two effects presented in section |5.2.1| due to the neutron 
distribution and the position of the reaction. Here only results of the study of the reaction d V 
are presented. The results of the target studies for the d +'^^ CI reaction are very similar. In order 
to study the effect of different neutron emission angles, two events with relatively high and low 
residual nucleus (^^Cr) kinetic energies were chosen from the evapOR simulation. The energies 
were i?RNi = 82.4 MeV and £'rn2 = 73.6 MeV in the laboratory system. The difference is due 
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to the neutron emission angle: the event with energy E^rni had a large neutron emission angle, 



— 174.4° while the event with £^rn2 had a small angle, 



^n,CM 



= 6.7°. The energies E] 



RNl 



and £'rn2 were used as input to TRIM. The energy losses of the residual nuclei are shown in figure 
[671 
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Figure 6.7: TRIM calculations of energy straggling of ^^Cr ions in a 200 /im/cm thick titanium 
target. In the figure, the kinetic energy distributions of the residual nuclei when they exit the target 
material are plotted for the case when neutrons are emitted at the angles 9^ — 174.4° (right) and 
CM = (left) corresponds to initial recoil energy of E'rni = 82.4 MeV and £'rn2 = 73.6 MeV, 
respectively. 



For the two cases the average energy losses are 

E^°^\ = 82.40 MeV - 78.01 MeV = 4.390 MeV, 



which gives the velocity 

/\„,loss 

c 

And 

which gives the velocity 



V rnjiMiC^ ^ ' 



1.5760 • 10" 



(6.1) 



TTlloSS 



73.60 MeV - 69.17 MeV = 4.430 MeV, 



A loss 



RN,in 



RN,out 



1.6858 • 10" 



(6.2) 



Next, the effects of this energy-loss difference on the Doppler shifts were studied. In the case 
of the ^^Cr residual nucleus, the 7-ray energy from the 2+ excited state to the 0+ ground state is 
1434 keV. For detector position at 90°relative to the incoming beam and a distance of 15 cm from 



the target to the front face of the detector, which was chosen for these simulations (see section 6.4 1 
the smallest angle for 7-ray detection, which corresponds to the largest Doppler shift, is about 70° 
Thus the difference in 7 ray energy becomes 



A loss Ai)^°^^ 



^^cosd^ ^ 0.054keV, 



(6.3) 



which is a small effect. The contribution from the Doppler shift from this studied effect is hence 
negligible. 

The second target effect concerns the position of the reaction. A new TRIM calculation was 
performed for ^^V ions with an energy oi E ^ 81.6 MeV. The energy loss of the ^^V beam in the 
target is shown in figure |6. 
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Figure 6.8: Energy distribution of ^^V residual nuclei for the case that the reaction occurs at the 
end of the target. The energy of the incident ^^V ion beam is 81.6 MeV. The target was a 200 
/xg/cm^ thick Ti foil. 

The energy loss of the beam through the target is, 

ElTy = 81.6 - 77.3 = 4.3 MeV. 

From the previous study of ^■^Cr, the average of -Erni and Eyw2 gives the average energy loss in 
the case the reactions occur at the beginning of target: 

. = EnNi+EnN2 ^ 4 ^^^^^ 
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The velocity of the ° Cr residual nucleus for the case of a reaction at the end of the target is 



- = 0.0571, (6.4) 

c 

and 

- 0.0570, (6.5) 

c 

if the reaction occurs in the beginning of the target. Note that it is assumed that the 7 rays are 
emitted after the residual nucleus has exited the target (see section ). The difference in velocities 
is small and gives a negligible effect on the Doppler shifts at 0^ ~ 70°: 

AE^ ^ E^„-^ cos 0^ = 0.05keV. (6.6) 

Hence, the position of the reaction in the target can be neglected. 

A similar study was done for the d('^''Cl, n) reaction showing a negligible effect on the Doppler 
shifts. The average energy loss for the d('^^Cl, n) reaction was 

ElTci = 3.50 MeV. (6.7) 

6.2.2 Angular straggling 



In figure 6.9 the angular deviation from the beam direction of the ^^V ions after passing through 
the 200 /ig/cm^ thick Ti target is shown. 
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Figure 6.9: Angular straggling of 81.6 MeV ^^V ions in a 200 jig/cm^ thick Ti target. 
The angular distribution of the ions due to straggling has a maximum at 0.1°. This is much 



smaller than the maximum of the angular distribution due to the neutron emission, see fig. 6.14 
Thus, the angular straggling in the target can be neglected. 



6.2.3 Life-time effects 

The compound nucleus has a very short life time compared to the time it takes for it to pass 
through the target. To prevent the 7 rays to be emitted while the residual nucleus is still in the 
target, it is required that the life time are long enough (see section |4|. If the 7 rays are emitted 
while the residual nucleus still is located inside the target a too large spread of Doppler shift will 
be the result. The average flight distances of the residual nuclei before emitting the 7 rays were 
studied. For the reaction d(^^V,n)^^Cr the average flight distance x, before the emission of the 
2+ —> 0+ 7-ray in ^^Cr is 

X « 1.0099ps • 0.0548c = 1.15 • lO^V = 16.5/im, (6.8) 

which is an order of magnitude larger than the target thickness w O.bfim. The life-time of the 
2+ 0+ transition in '^^Ar is 0.578 ps and gives 



X « 0.578ps • 0.0589c = lO.O^im, (6.9) 

which is also large compared to the target thickness. Thus the photons will be emitted outside the 
target material in both reactions. 

The effective life time of the transitions of interest must also be shorter than the time it takes 
for the residual nuclei to travel a distance < 1mm. This effect was not studied in the present work. 



6.2.4 Experimental determination of the intrinsic FWHM 



The intrinsic part of the FWHM, Wi, given in equation 3.10 is unique for every detector type. 
This value is used as input to the tracking program MGT. Data from a measurement of the core 
segment of AGATA crystal number C002 using a -^^^Ra source |34] was used to obtain the energy 
dependence of Wi for a typical AGATA detector. A straight line was fitted to the data points, see 
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Figure 6.10: Intrinsic FWHM, Wi, as a function of 7-ray energy for AGATA crystal C002 "M]. 



Using the fitted straight line of fig. |6.10| the intrinsic FWHM of the peaks corresponding 
to the 1434keV and 2168keV 2+ 0+ transitions in ^^Cr and ^^Ar are 2.30keV and 2.83keV, 
respectively. 



6.3 evapOR 

The studies of the target effect described in section |6.2[ yielded an average energy loss in the 
target of 4.3MeV and 3.5MeV for the ^^V and '^^Cl induced reactions, respectively. The laboratory 
energies selected for the d(^^V,n)^^Cr and d('^^Cl, n)'^*Ar reactions were 77.0 Me V and 66.5 MeV, 
respectively. 

6.3.1 Energy and angular distributions 

Figure |6.11| shows the energy and angular distribution of the neutrons as obtained by evapOR. 
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Figure 6.11: The left panel shows the center of mass neutron energy and the right panel 

shows the angular distribution of the neutrons in the laboratory system, 0£ab '^or the reaction 
d (Sly, n) Cr at E (^i V) = 77MeV. 



The energy distribution (left hand figure) is Maxwellian as expected from equation 2.11 The 
distribution has a maximum at about IMeV. According to the plot on the right hand side of 
figure |6.11| the neutrons have a large emission angle in the laboratory frame. The centroid of the 
distribution is around 40°. 

Figure |6.12| shows the angular distribution of the neutrons in the center of mass system. 



Neutron Angle distribution in CM frame 
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Figure 6.12: Angular distribution of neutrons in the center of mass system for the reaction 
d(^"'^V, n)^^Cr at E (^"'^V) = 77MeV. The black histogram shows the number of neutrons emitted 
per OcM value. The red histogram is the black histogram divided by s'm {6cm) .The histograms 
are arbitrarily normalized to each other. 
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Figure 6.13: Velocity distributions of the residual nuclei for the d( Cl,nj Ar (left) and 



d (^iV,n)^^ Cr (right) reactions at E (^^Cl) = 66.5MeV and E (^ly) = 77MeV, respectively. 



Figure 6.13 shows the velocity distribution of the residual nuclei in both reactions. 



The neu- 
The 



tron energy distributions in the laboratory frame are shown in the left panel of figure 6.14 
d(^"'^V, n)^^Cr reaction contains more neutron events and thus the area under the red histogram is 
larger. The right panel of figure [6?T4| shows the angular distributions of the residual nuclei. There 
is an evident shift of the angular distribution, implying a smaller mean deflection of the ^^Cr nuclei 
compared to the '^^ Ar nuclei. This difference in angular deflection plays a major role in the Doppler 
shift of the 7-ray energies. 




Figure 6.14: The left panel shows the energy distributions of the neutrons and the right panel the 
angular distribution of the residual nuclei in the laboratory system. 



6.4 Simulated 7-ray spectra 

6.4.1 Graphical illustration of the GEANT4 simulation 



In figure 6.15| a graphical illustration of the Geant4 simulation is shown. Ten events of the 



reaction ^^V (77 MeV) + d were generated by evapOR. The direction of motion of the incoming 
^^V beam is from right to left in the figure. For the plot only events belonging to the reaction 
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Figure 6.15: A graphical illustration of a Geant4 simulation in which a compound nucleus emits 
neutrons and 7 rays, which are impinging on an AGATA triple cluster detector. See text for further 
details. 



The 7 rays were forced to be emitted isotropically into a cone with an opening angle of 20°, 
centered at 90° relative to direction of incoming beam. An AGATA triple cluster detector was 
placed at that angle, and covered a solid angle which was smaller than the solid angle of the 7-ray 
emission cone. The distance between the 7-ray source and the front of the HPGe crystals was 15 
cm in the figure. 

Except if otherwise noted, the setup shown in fig. 6.15| with the angle and distance given 
above, was used in the simulations presented in the following subsections. 



6.4.2 Gamma-ray spectrum produced by summing 

As a first simple test, a 7-ray spectrum was created by summing for each event the total energy 
deposited in the detector. Such a spectrum is shown in fig. |6.16| for the 2168 keV 7 ray emitted 
in the d('^''Cl, n)^^Ar reaction at i?LAB ('^^Cl) — 66.5 MeV. The 7 rays were emitted isotropically 
into a 90° ± 0.5° cone. The intrinsic energy resolution was set to {Wi — keV) and no Doppler 
corrections were applied. The full-energy peak, the single- and double-escape peaks, and Compton 
background are clearly visible in the spectrum. 
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Figure 6.16: A 7-ray spectrum obtained by summing the total energy absorbed in the AGATA 
triple cluster detector. A total of 10^ evapOR events of the reaction '^^Cl(66.5 MeV) + d were 
generated. See text for further details. 



6.4.3 Gamma-ray spectra produced by tracking and Doppler correction 

Tracked 7-ray energy spectra were created by feeding the output of the Geant4 simulation into 
the MGT tracking program. The MGT smearing parameter, which emulates a finite interaction 
position resolution in the HPGe crystals, was used as one of the parameters in the simulations. 

The 7-ray emission angle, which is needed by the Doppler correction procedure, was calculated 
by assuming that the 7 rays were emitted from the center of AGATA and detected at the first 
interaction point, which was obtained from MGT. For the Doppler correction the velocity vector of 
the residual nucleus is also needed. Except were otherwise noted, an average value of the velocity 
vector was used, namely an angle of 0° (motion parallel to the beam) and the average velocity 
given by the evapOR simulation. The use of such average values is common in fusion-evaporation 
reactions, when no ancillary detector is available for the event-by-event detection of the velocity 
vector of the residual nuclei. 

Initially, two simulations were performed to verify that the tracking and Doppler correction 



procedures functioned as expected. In the first simulation, shown in the left panel of fig. 6.17 the 
1434 keV 7 rays from the reaction d (^^V, n)^^ Cr, at i?LAB (^^V) = 77.0 MeV, were emitted at an 
angle of exactly 90° relative to the incoming beam. Gamma-ray tracking with a 5 mm smearing 
was used, but no Doppler corrections were applied. In the second simulation, shown in the right 
panel of fig. |6.17| the 1434 keV 7 rays were emitted isotropically into the 90° ± 20° cone covering 
the solid angle subtended by the triple cluster detector, and both tracking (with 5 mm smearing) 
and Doppler corrections were applied. The distance from the source to the front of the HPGe 
crystals was 15.0 cm. The FWHM of the 1434 keV peak in both panels of fig. 6.17 are almost 
identical, which shows that the 7-ray tracking and Doppler correction algorithms functioned as 
expected. 
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Reaction 


E^[keV] 


W [keV] 


Wi [keV] 


Wr, [keV] 


d(37Cl,n)'"Ar 


2168 


7.74 


2.83 


7.20 


d(5iV,n)'' Cr 


1434 


4.55 


2.30 


3.93 



Table 6.2: Contributions to the total FWHM, W, due to the intrinsic resolution, Wi, and the 
Doppler effects, Wb, for the two studied reactions. 




Figure 6.17: Spectra of tracked 7 rays emitted in the reaction d(^^V, n)^^Cr at E'lab (^^V) = 77.0 
MeV. Left: The 7 rays were emitted at exactly 90° relative to the incoming beam and no Doppler 
correction was applied. Right: The 7 rays were emitted isotropically into a 90° ± 20°cone covering 
the solid angle subtended by the triple cluster detector, and a Doppler correction was applied. See 
text for further details. 



Results of tracked and Doppler corrected 7-ray spectra from simulations in which also the 



energy dependence of the intrinsic energy resolution (see section 6.2.41 has been included, are 



shown for the two reaction in fig. 6.18 and 6.19 The smearing parameter was kept at 5 mm and 
the 7 rays were emitted isotropically into the 90° ± 20° cone. 
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Figure 6.18: Tracked and Doppler corrected 7-ray spectrum of the 1434 keV 7 rays emitted in the 
reaction d ^^^Vju)^^ Cr. See text for further details. 
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Figure 6.19: Tracked and Doppler corrected 7-ray spectra of the 1434 keV (left) and 2168 keV 
(right) 7 rays emitted in the reactions d{^'^Y ,n)'^'^Cr {Ei^ab (^^V) = 77 MeV) and d(3^Cl, n)38Ar 
{EhAB ("^^Cl) — 66.5 MeV), respectively. See text for further details. 



Numerical values of the total FWHM W of the peaks in fig. 6.19 and their components (see 



eq. 3.9 1 due to the intrinsic resolution {Wi obtained from fig. 6.10) and due to the Doppler eff'ects 
Wb, are given in table |6.1 



The results given in the table clearly show that the contribution to the total FWHM, due to the 
Doppler efi^ects, is much larger in the reaction d ('^^Cl, n)^^ Ar than in the reaction d (^^V, n)^^ Cr. 
The main reason for this is the difference in mass number of the compound nuclei '^^Ar and ^^Cr. 
This mass difference leads to a somewhat larger and broader velocity distribution of the residual 
nuclei (see fig 6.131 and, more importantly, to an angular distribution of the residual nuclei which 
is both broader and has a maximum at a larger angle (see right panel of fig. 6.14) for the reaction 
d (^^Cl,n) Ar compared to d (^""^Vju) Cr. The angular and velocity distributions of the 7-ray 
emitting residual nuclei infiuence directly the Doppler effects, which is seen as an effect on the 
FWHM of the peaks. 

A simulation was also performed by correcting for the Doppler effects introduced by the varia- 
tion the velocity vector of the residual nucleus. A comparison of the FWHM of the 1434 keV peak 
obtained when using the average velocity v ector and the precise value, as obtained event-by-event 

For the reaction d (^^V,n) at Slab (^^V) = 77.0 



6.20 



from the evapOR simulation, is shown in fig 
MeV the FWHM is improved from 4.7 keV to 3.9 keV when the precise value of the velocity vector 
of the residual nucleus is applied in each event. 
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Figure 6.20: Comparison of the FWHM of the 1434 keV peak when the Doppler correction was 
performed by using for the residual nucleus an average velocity of v/c = 0.0548 and an angle 
of 0° (left) and the exact value obtained event-by-event from evapOR (right). The reaction was 
d(^"'^V,n) Cr at i?LAB (^^V) — 77.0 MeV, the smearing parameter was 5 mm. 



6.4.4 Doppler effects induced by emission of a particles 

A simulation of Doppler effects, following a reaction channel in which a particles are emitted, is 
shown in fig. |6.21[ The emitted a particles are four times heavier than neutrons and have much 



higher kinetic energies in the laboratory system (compare the left panels of fig. 6.14 and 6.211, 
which leads to very large spreads both of the energy and angle of the residual nuclei. This in turn 
leads to very broad 7-ray peaks, as seen in the right panel of fig. |6.21[ Thus, for reaction channels 
with emission of a particles it is usually necessary to determine event-by-event the velocity vector 
of the residual nuclei. 
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Figure 6.21: Doppler effects following the reaction d(3^Cl,a)35S at £:lab(^^C1) = 66.5 MeV. Left: 
Kinetic energy distribution of the residual nuclei ^^S. Right: FWHM of a 7-ray peak at 2168 keV 
after tracking using the average residual nucleus velocity vector for the Doppler correction . The 
smearing parameter was 5 mm. 



6.4.5 FWHM as a function of smearing parameter 

Simulations, in which the total FWHM of the 1434 keV and 2168 keV peaks were determined as 
a function of the MGT smearing parameter, are shown in figures [6.22| and |6.23| for two different 
detector distances, 15.0 cm and 23.5 cm, respectively. As seen in the figures, the slopes increase for 
small values of the smearing parameter, up to about 4-8 mm and become more or less constant for 
larger values. The slopes are quite similar in all cases, although slightly larger for the d(^^ V, n)^^Cr 
reaction and for the shorter distance of 15.0 cm. 
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Figure 6.22: FWHM as a function of the MGT smearing parameter for the 2168 keV (left) and 
1434 keV (right) peaks of the d(^'''Cl, n)^*Ar and d^^^ V, n)^^Cr reactions, respectively. The distance 
from the source to the front of the HPGe crystals was 15.0 cm. The data points are connected by 
straight lines. 




Position resolution (mm) Energy (keV) 

Figure 6.23: Same as fig. |6.22| but for a distance of 23.5 cm from the source to the front of the 
HPGe crystals. 



The interaction position resolution which is expected for the AGATA HPGe crystals is of the 
order of 5 mm, a value which, however, is not constant. The resolution will depend on the position 
of the interaction in the crystal. For example in the front part of the crystal, where the electric 
field is rather non-uniform, the interaction position resolution will be worse compared to a position 
in the central parts. The interaction position resolution will also depend on the energy deposited 
in the interaction point. It is expected that it will be proportional to the inverse square root of the 
interaction energy, a dependence which is implemented in MGT, although not used in this work. 

A conclusion of the results shown in fig. |6.22| and |6.23| is that with the proposed reactions, 
d(37Cl,n)38Ar at £;lab(^^C1) = 66.5 MeV and d(5iV,n)52Cr at £;lab(^^V) = 77.0 MeV, it will 
be difficult to get a good determination of the interaction position resolution from the measured 
FWHM of the peaks. One possibility to improve the sensitive, is to make the measurements 
for much smaller distances, than what was used in this work (15.0 cm and 23.5 cm). Another 
possibility is to try to find another reaction with larger average recoil velocities but which still has 
small enough spread of the residual nucleus velocity vector so that the average Doppler correction 
procedure may be used. 
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6.5 Conclusions and summary 



In section |6.1| the studied reactions in TALYS resulted in two interesting reactions for further 
study, namely d(^^V, n)^^Cr and d('^^Cl, n)'^*Ar. On the basis of the requirements in section|4] these 
reactions have suitable 7-ray energies and life times within the required range. Both reactions have 
dominant cross sections for the desired neutron channel. Both reactions lead to even-even residual 
nuclei and give quite large recoil velocities (important in order to study the Doppler effects). 

In section [6^ the target effects were investigated using TRIM. Two major sources of unwanted 
target effects were studied and determined to be neghgible for the Doppler shifts. The average 
energy loss for the ions in the target material were simulated and presented. 

In section 6.3 results from evapOR and TALYS were compared. Different energy and angular 
distributions were produced by evapOR. Due to the different amount of nucleons in the compound 
nuclei of the two reactions, the ^®Ar residual nuclei were deflected more from the ion beam direction 
compared to ^^Cr. This difference in angular deflection is the reason why the Doppler broadening 
is larger for the reaction d('^^Cl, n)'^^Ar than for d(^^V, n)^^Cr. 

In section 6.4 the detection of the 7 rays in one AG ATA triple cluster detectors was simulated 
with Geant4. The 7 rays were tracked by the MGT tracking program to produce Doppler cor- 
rected 7-ray spectra. The FWHM of the 7 ray peaks of interest were calculated as a function of the 
interaction position resolution. The d(^^V,n)^^Cr reaction was found to be slightly more suitable 
than the d (^''Cl, n) Ar reaction for the planned AGATA commissioning experiment. 
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Appendix B 



Available ion beams at Laboratori 
Nazionali di Legnaro (Padova) 



Table B.l: Available ion beams at LNL 



Beam 


Current in (nA) 


Beam 


Current in (nA) 


Beam 


Current in (nA) 


IH 


1000 


40Ca 


150 


79Br 


800 


2H 


300 


48Ca 


150 


81Br 


800 


6Li 


50 


48Ti 


700 


90Zr 


100 


7Li 


200 


50Cr 


150 


91Zr 


25 


lOB 


100 


52Cr 


150 


92Zr 


30 


UB 


400 


51V 


250 


94Zr 


30 


12C 


3000 


54Fe 


200 


96Zr 


300 


13C 


35 


56Fe 


200 


92Mo 


250 


14N 


800 


58Ni 


1500 


94Mo 


150 


160 


2000 


60Ni 


500 


95Mo 


250 


170 


500 


64Ni 


500 


96Mo 


250 


180 


500 


63Cu 


1000 


97Mo 


150 


19F 


3000 


65Cu 


350 


98Mo 


400 


24Mg 


300 


64Zn 


500 


lOOMo 


150 


26Mg 


200 


66Zn 


250 


96Ru 


200 


27A1 


400 


68Zn 


200 


98Ru 


70 


28Si 


1000 


69Ga 


300 


99Ru 


450 


29Si 


300 


70Zn 


200 


lOORu 


450 


30Si 


200 


71Ga 


200 


lOlRu 


600 


31P 


500 


74Ge 


800 


102Ru 


1000 


32S 


2000 


76Ge 


200 


104Ru 


650 


33S 


400 


76Se 


200 


107Ag 


400 


34S 


200 


77Se 


300 


109Ag 


400 


36S 


400 


78Se 


1000 


1271 


800 


35C1 


2000 


80Se 


2000 


197Au 


500 


37C1 


650 


82Se 


400 
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Appendix C 

Simulation steps 



1. evapOR 

• Run cvapOR. 

• Modify the *.inp file for evapOR, where you can choose nuclei, spin and energy. 

• evapOR-loop 

• evapOR-loop *.inp 

You will have now some produced files. The important one is *.pax 

2. Sortpax: 

You run sortpax to generate the Doppler shifted gammas and to produce the *.aga file. 

• sortpax *.pax 

• choose aga file with typing '*.aga' 

• choose txt file to read the events, with typing '*.txt' 

• choose bin file for evapOR Root tree, with typing '*.bin' 

• choose the gate on nucleus (Z, N) 

• choose the angle of emitted photons (thetamin, thetamax, phimin, phimax) 
Now you will have the *.aga file. (This you copy to fari in my case) 

3. AGATA Geant4: 

Go to directory of AGATA Geant4 

• Run 'Agata -Ext -n' in order to read the particles also. 

• Type the desired of the following commands 

— /Agata/file/enableLM 

— /Agata/file/verbose 1 to write information 

— /Agata/file/packingDistance 5. for packing 

• Define the geometry, a cluster is given by this macro 

— /control/execute macros/geomlSOAli.mac 
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— Look in macro file and manual agata Geant4 to see the other geometries. The distance 
to the detector is given in AlSOeulerAli.list which is under the folder A180. You can 
there put and move the detector anywhere. 

To generate with the *.aga file use 

• /Agata/generator/emitter/eventFile /homc/tslO/ali/Desktop/dopplcr.aga /Agata/run/beamOn 
100000 

To generate with AGATA 

• /Agata/generator/recoil/beta 5. (for nucleu's velocity) 

• /Agata/generator/gamma/energy 1000 (the energy of gammas) 

• /Agata/gcncrator/gamma/thctaRangc 65.0 115.0 (the angle emitted gammas) 

• /Agata/generator/gamma/phiRange 65.0 115.0 (the </> angle emitted gammas) 

4. MGT: 

The output is gammaevents. 0000/1/2 and this should be used in MGT go the folder where this 

file is and type 

• ' mgt -fw -f GammaEvents.OOOO -11 -oa 30 30 -dd 5 5 -sr 0.05 ' 

where fw is defining the intrinsic FWHM, statistical and then noise factor. 11 0/1/2/3 gives more 
information, dd defines the smearing and packing . 

• sr gives the average velocity and direction in 9 and (f). 
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